Polymorphic Solid Solutions in Molecular Crystals: Tips, Tricks, and Switches.

Crystal polymorphism has been a topic of much interest for the past 20 years or so, especially since its scientific (and legal) importance to the pharmaceutical industry was realized. By contrast, the formation of solid solutions in molecular crystals has been overlooked despite its long-standing prevalence in the analogous field of inorganic crystals. Wilfully forgotten, crystalline molecular solid solutions may be very common in our world since molecular compounds are rarely produced with 100% purity, and impurities able to form solid solutions are difficult to reject via recrystallization. Given the importance of both polymorphism and solid solutions in molecular crystals, we share here some tips, tricks, and observations to aid in their understanding. First, we propose a nomenclature system fit for the description of molecular crystalline solid solutions capable of polymorphism (tips). Second, we highlight the challenges associated with their experimental and computational characterization (tricks). Third, we show that our recently reported observation that polymorph stabilities can change by virtue of solid solution formation is a general phenomenon, reporting it on a second system (switches). Our work focuses on the historically important compound benzamide forming solid solutions with nicotinamide and 3-fluorobenzamide.

Rotamers in Crystal Structures of Xylitol, D-Arabitol and L-Arabitol.

Rotamers are stereoisomers produced by rotation (twisting) about σ bonds and are often rapidly interconverting at room temperature. Xylitol-massively produced sweetener-(2R,3r,4S)-pentane-1,2,3,4,5-pentol) forms rotamers from the linear conformer by rotation of a xylitol fragment around the C2-C3 bond (rotamer 1) or the C3-C4 bond (rotamer 2). The rotamers form two distinguishable structures. Small differences in geometry of rotamers of the main carbon chain were confirmed by theoretical calculations; however, they were beyond the capabilities of the X-ray powder diffraction technique due to the almost identical unit cell parameters. In the case of rotamers of similar compounds, the rotations occurred mostly within hydroxyl groups likewise rotations in L-arabitol and D-arabitol, which are discussed in this work. Our results, supported by theoretical calculations, showed that energetic differences are slightly higher for rotamers with rotations within hydroxyl groups instead of a carbon chain.

First Experimental Quantitative Charge Density Studies of Advanced Intermediate of Vitamin D Analogues.

Vitamins D are a group of fat-soluble secosteroids which play a regulatory role in the functioning of most cells. Rational design of new vitamin D analogs, of increased therapeutic potency and lowered calcemic side effects, requires high-resolution initial structures and a deep understanding of interactions with the molecular targets. In this paper, using quantum crystallography, we present the first determination of the experimental quantitative charge density of an advanced intermediate of vitamin D analogues as well as a reconstruction of the theoretical electron density of final vitamin D analogues. Application of these methods allows for topological and electrostatic interaction energy analysis. We showed that the A-ring chair conformation has a significant influence on the topological properties of vitamin D compounds. Moreover, the interactions between the CD-ring and side-chain additionally stabilize the crystal structure. These results are supported by our theoretical calculations and previous biological studies.

Further Validation of Quantum Crystallography Approaches.

Quantum crystallography is a fast-developing multidisciplinary area of crystallography. In this work, we analyse the influence of different charge density models (i.e., the multipole model (MM), Hirshfeld atom refinement (HAR), and the transferable aspherical atom model (TAAM)), modelling of the thermal motion of hydrogen atoms (anisotropic, isotropic, and with the aid of SHADE or NoMoRe), and the type of radiation used (Mo Kα and Cu Kα) on the final results. To achieve this aim, we performed a series of refinements against X-ray diffraction data for three model compounds and compared their final structures, geometries, shapes of ADPs, and charge density distributions. Our results were also supported by theoretical calculations that enabled comparisons of the lattice energies of these structures. It appears that geometrical parameters are better described (closer to the neutron values) when HAR is used; however, bonds to H atoms more closely match neutron values after MM or TAAM refinement. Our analysis shows the superiority of the NoMoRe method in the description of H-atom ADPs. Moreover, the shapes of the ADPs of H atoms, as well as their electron density distributions, were better described with low-resolution Cu Kα data in comparison to low-resolution Mo Kα data.

Relation Between Crystal Structures of Precursors and Final Products: Example of Vitamin D Intermediates.

In this paper, we proved that the solid-state structure of vitamin D analog is well represented by the structures of its structural fragments. This is important in predicting the biological activity of vitamin D analogs that are not available in the solid form. The previously published crystal structure of advanced vitamin D intermediate provided additional insights into vitamin D properties. A similar analysis based on simple vitamin D intermediate analogues showed that precursors crystallized in the space groups typical for vitamins D; geometrical parameters were related to the corresponding parameters in the vitamin D analogues; and crystal structures of the basic intermediates and their final products contained similar intermolecular interactions, essential for the infinite hydrogen bond motif observed in the vitamin D analogues. The energy of these interactions is related as shown by theoretical calculations, that is, energy frameworks analysis. Moreover, analysis of the hydrogen bonds motifs revealed a relation between these motifs and the absolute configuration of basic intermediates as well as the space orientation of the exocyclic methylene group in the final structures.

Energetics of interactions in the solid state of 2-hydroxy-8-X-quinoline derivatives (X = Cl, Br, I, S-Ph): comparison of Hirshfeld atom, X-ray wavefunction and multipole refinements.

In this work, two methods of high-resolution X-ray data refinement: multipole refinement (MM) and Hirshfeld atom refinement (HAR) - together with X-ray wavefunction refinement (XWR) - are applied to investigate the refinement of positions and anisotropic thermal motion of hydrogen atoms, experiment-based reconstruction of electron density, refinement of anharmonic thermal vibrations, as well as the effects of excluding the weakest reflections in the refinement. The study is based on X-ray data sets of varying quality collected for the crystals of four quinoline derivatives with Cl, Br, I atoms and the -S-Ph group as substituents. Energetic investigations are performed, comprising the calculation of the energy of intermolecular interactions, cohesive and geometrical relaxation energy. The results obtained for experimentally derived structures are verified against the values calculated for structures optimized using dispersion-corrected periodic density functional theory. For the high-quality data sets (the Cl and -S-Ph compounds), both MM and XWR could be successfully used to refine the atomic displacement parameters and the positions of hydrogen atoms; however, the bond lengths obtained with XWR were more precise and closer to the theoretical values. In the application to the more challenging data sets (the Br and I compounds), only XWR enabled free refinement of hydrogen atom geometrical parameters, nevertheless, the results clearly showed poor data quality. For both refinement methods, the energy values (intermolecular interactions, cohesive and relaxation) calculated for the experimental structures were in similar agreement with the values associated with the optimized structures - the most significant divergences were observed when experimental geometries were biased by poor data quality. XWR was found to be more robust in avoiding incorrect distortions of the reconstructed electron density as a result of data quality issues. Based on the problem of anharmonic thermal motion refinement, this study reveals that for the most correct interpretation of the obtained results, it is necessary to use the complete data set, including the weak reflections in order to draw conclusions.